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The last two years have seen an impressive improvement in the determination of \V u b\, especially from inclusive 
decays. The error on \V u b\ measured with inclusive decays was reduced from 18% (PDG 2004) to 8% (PDG 2006). 
This progress is a result of combined experimental and theoretical efforts. In this talk, the theoretical framework 
(BLNP) that enabled such progress is reviewed, as well as other approaches to an inclusive determination of \V u b\ 
(DGE, Mx — q 2 etc.). The prospects of improving \V u t\ are discussed, addressing issues of weak annihilation, 
implications of leptonic B decays, and determination of \V u b\ with exclusive decays. 



1. INTRODUCTION 

V u b i one of the smallest matrix element of the 
CKM matrix, is one of the fundamental param- 
eters of the Standard Model. In the geometrical 
picture of the unitarity triangle, the side opposite 
to the angle (3 (= (f>\ ) is proportional to | V u b/V c b \. 
Since this angle is measured to high accuracy by 
the B factories, and the error on \V c b\ is at a level 
of 2% [1], an accurate measurement of \V u b\ is 
important for constraining the unitarity triangle. 

\V u b\ can be measured either through exclu- 
sive decays, e.g. B — > ttI~9, or inclusive decays: 
B — > X u l~v. In the former method we encounter 
a large theoretical uncertainty as a result of our 
limited knowledge of the form factors that control 
the decay. The latter method offers, in principle, 
the most accurate way to extract \V u b\ from the 
total width r (B — > X u l~D). In practice, since 
I Kb | \Vub\, the total width cannot be mea- 
sured due to the large charm background from 
B — > X c l~D decays. In order to eliminate the 
charm background one is forced to look at re- 
gions of phase space where charm cannot be pro- 
duced. In these regions the theoretical descrip- 
tion is more complicated, but thanks to recent 
progress in our understanding of charmless inclu- 
sive B decays, the inclusive measurement of \V u b\ 
gives a smaller error compared to the exclusive 
one. 

The recent global analysis of the unitarity tri- 
angle shows a "tension" at the level of 2a between 



(= 4>i) and \V u b\ [2] (see the contribution by S. 
T'Jampens at this conference). This is a result 
of two effects. One, the central value of \V u b\ has 
increased, especially as measured from inclusive 
decays, and the error bar decreased from 18% in 
2004 [3] to 8% in 2006 [1]. At the same time 
the value of sin2/3(= sin20i) has decreased. It is 
therefore important to understand how the value 
of \V u b\ is obtained. 

In this talk we will therefore focus on the the- 
ory of inclusive measurement of \V u b\. For various 
experimental issues wc refer the reader to the con- 
tribution by L. Gibbons at this conference, and 
for exclusive measurement using lattice data we 
refer to the contribution of C. Davies. 

2. KINEMATICS 

We begin our discussion by shortly reviewing 
the kinematics of semileptonic B decays. For a 
"pedestrian" introduction to semileptonic B de- 
cays see chapter 1 of [4]. 

Any B — > X u l~v event can be described by 
three kinematical variables. The triple differ- 
ential decay rate depends, via the optical theo- 
rem, on the hadronic tensor (for derivation 
see [4]). Given a basis of two four vectors, the 
hadronic tensor is usually decomposed into the 
various possible Lorentz structures, where the co- 
efficients of these structures are called structure 
functions and denoted as W%. 

In [5] the use of the (v, n) basis was advocated. 
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In this basis v is the four velocity of the decay- 
ing B meson and n is a light-like vector in the 
direction of the hadronic jet. In the rest frame 
of the decaying meson we have v — (1,0,0,0) 
and we can choose n to be n = (1,0,0,1). One 
can also define a conjugate light like vector n — 
(1, 0, 0, —1), such that 2v — n + n (in the follow- 
ing we assume these values for v, n and h). This 
choice of vectors motivates the following choice of 
kinematical variables: 



P l = M B -2E i , 
P-=n-P x =E x + \P x \ 
P+=n-P x =E x -\P x \ 



(1) 



where P x is the four momentum of the hadronic 
jet and Ei is the energy of the charged lepton. 
Other kinematical variables can be expressed in 
terms of this choice of variables. For example 
the hadronic and leptonic invariant masses are: 
M\ = P+P_, and q 2 = (M B - P_)(M B - P+), 
respectively. 

In terms of these variables the exact expression 
for the triple differential decay rate is: 



d 3 T„ 



dP+ dP- dPi 
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(M B -P+) 



(P_ - Pi){M B - P- + Pi - P+) W 1 
+ (M B -P_)(P_-P + )^ + 
(P_ - P ; )(P - P+) ( |Ws + Wi + -W 5 

v y 



(2) 



where y = (P_ - P + )/(M B - P + ) and Wi are de- 
fined in [5] and do not depend on p. This choice 
of variables and basis has two main advantages. 
The first is that the phase space has probably the 
simplest form possible: 



Ml 



< P+ < Pi < P- 



< M f 



(3) 



The second advantage is that there is no explicit 
dependence on the mass of the b quark in the 
expression for the triple rate (2). This allows for 
theoretical predictions of partial rates instead of 
event fractions and eliminates the large source of 



uncertainty from the value nib- The triple rate 
depends on nib only through the W functions, 
which we now discuss. 

3. DYNAMICS 

The structure functions Wi cannot be calcu- 
lated exactly. Fortunately for heavy meson de- 
cays there are two small parameters we can ex- 
pand Wi in: the mass of the b quark (or more 
exactly Aqcd/ rnb) and a s . 

If we had no charm background we could inte- 
grate over P± up to M B and use a Heavy Quark 
Effective Theory (HQET) based^ Operator Prod- 
uct Expansion (OPE) to write W as: 



+ 



(4) 



(We use a short hand notation in here. In practice 
there are several operators at each order and the 
coefficients are generalized functions and not con- 
stants). The coefficients Ci(pi) contain the short 
distance physics (p, ~ nib) and are calculable 
in perturbation theory ((Oi) = as a result of 
HQET equations of motion [6] ) . Currently Co [7] 
is known at 0(a s ) while Ci [8,9] and C3 [10] are 
known at 0(a°). Recently even C4 was calculated 
at O(a s ) [11]. 

The matrix elements of the local operators Oi 
between B meson states, (Oi), are called the 
Heavy Quark (HQ) parameters. They contain the 
long distance physics and must be taken from ex- 
periment. We expect them to scale as (Oi) ~ 
Aqcd- We nave (^o) = 1; (O2) defines two HQ 
parameters n 2 and a 2 , = 3[(M B ) 2 - (M B ) 2 ]/A; 
(O3) defines p\ s and p z D etc. 

A similar OPE can be constructed for B — > 
X c l~v and B — ► X s j, which contains the same 
HQ parameters. For B — > X c l~v the OPE works 
very well. As a result the error on |V^&| is at the 
level of 2% [1] and the HQ parameters can be ex- 
tracted from experiment. For B — > X s j a local 
OPE is more problematic for two reasons. The 
first is that the OPE is valid for photon energies 
that cannot be attained at the current experi- 
ments. The second is that even for such ener- 
gies when one goes beyond leading order in I /nib 
non-local operators arise (see [12] and references 
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within). 

In practice we have to deal with a large charm 
background from B — > X c l~v. In order to elim- 
inate this background experimental cuts must be 
imposed. The "charmless" region is typically the 
region of phase space where P + ~ Aqcd and 
P- ~ mb ■ In this region, since P + -C P_, the 
HQET based OPE is not valid (a fact commonly 
referred to as "breakdown of the OPE"). This re- 
gion is known as the shape-function (SF) region, 
to distinguish it from the region of P + <~ P_ ~ 
mb, the OPE region. 

Still, we do have a systematic 1/rrib expansion 
in this region. The Wi can be written as: 

Wi~H u - J®S+— Y,h k u -j k u ®s k u + ---. (5) 
mb k 

The H, J, and 6"s are called the hard, jet, and 
shape functions respectively, and they encode the 
physics at the scales fi^ ~ m,b, fa ~ y/rrib Aqcd, 
and no ~ Aqcd, respectively. The structure 
functions are sensitive to the "intermediate" scale 
fa ~ \/ m bJ^QCD since this was introduced as a re- 
sult of the experimental cuts. The hard and jet 
functions are calculable in perturbation theory, 
while the shape functions are non-perturbative 
objects. Currently the leading order hard func- 
tion is known at 0(a s ) [13,14] and the leading 
order jet function at order 0{a 2 s ) [15]. The sub- 
leading functions are only known at 0(a®) [16- 
18]. 

For the photon spectrum in B — > X s j near the 
kinematical end-point, which is the region mea- 
sured well by experiments, we have a similar ex- 
pansion: 

k 

Notice that the leading order jet and shape func- 
tion are the same as B — ► X u while the hard 
function is different. Beyond leading order this is 
no longer the case. 

What is the relation between the two regions? 
Neglecting a s corrections, moments of the shape 
functions are related to the HQ parameters. For 
example, the first moment of S, the leading or- 
der shape function, is related to mb, while the 



second moment is related to ^ [19]. Including 
a s corrections, we can define a renormalization 
scheme, called the "shape function scheme" [14], 
such that these relations hold at each order in 
perturbation theory. Currently, these relations 
are known at order a 2 s for mb and /j, 2 . [20]. The 
high precision of these relations implies that a 
good knowledge of the HQ parameters helps to 
constrain the shape functions. There are similar 
relations between the subleading shape functions 
and the HQ parameters, although they are known 
only at order a®. Still, such relations help us to 
model the subleading shape functions. 

This concludes our brief review of the dynamics 
of inclusive B decays. We should emphasize that 
the above description is not a theoretical model, 
but a rigorous theory based on QCD and a sys- 
tematic expansion in 1/rnb- Any inclusive extrac- 
tion of \V u b\ has to be based on these ingredients. 

4. INCLUSIVE EXTRACTION OF V ub - 
PRESENT 

There are currently several theoretical cal- 
culations, already implemented by experiments, 
which combine some or all of the ingredients de- 
scribed in the previous section. We now review 
each of them briefly (the LLR approach will be 
discussed in the next section). For a more de- 
tailed account see the original papers. 

BLNP Approach: The BLNP (Bosch-Lange- 
Neubert-Paz) approach [5] is a culmination of a 
research efforts that extends over 12 years. Its 
"philosophy" is to use all that is currently known 
about the triple differential decay rate of B — ► 
X u l~v and B — > X s ~/, namely 

• At leading order in l/m*,: H u , H s , J at 
0(a„); 

• 1/rrib subleading shape functions at O{a s ); 

• Known a s /rrib terms from the OPE calcu- 
lation (the part of Co in equation (4) that 
becomes subleading in the SF region); 

• Known \jm\ terms from the OPE calcula- 
tion (part of c 2 in equation (4)). 

The various ingredients are implemented in such a 
way that there is a smooth transition between the 
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SF and OPE regions, in the sense that once the 
kinematical variables are integrated far enough 
the OPE result is recovered. 

Non-perturbative physics effects are contained 
in the leading and subleading order shape func- 
tions. BLNP uses the fact that the leading or- 
der shape function is universal, so the formalism 
allows the extraction of the leading order shape 
function from the photon spectrum in B — > X s j 
and its use as an input for B — ► X u l~D (be- 
yond leading order a slight redefinition of of S 
is needed, see [5]). The subleading shape func- 
tions arc modeled in accordance with their mo- 
ment constraints. 

The complete error analysis consists of the fol- 
lowing: 

• The leading order shape function is taken 
from experiment as explained above, sub- 
stantially reducing the associated error. 

• The perturbative error is estimated by vary- 
ing the scale of a s in the different terms. 

• The subleading shape function error is esti- 
mated as follows. At order there are 3 
subleading functions and 9 models are con- 
structed for each of them. The 9 3 = 729 
possible combinations are then scanned to 
estimate the resulting error. The models 
for the subleading shape functions are con- 
structed such that they respect the moment 
constraints (the zeroth moment vanishes, 
the first moment is related to fj, 2 and p? G , 
and the second moment is of order Aq CD ), 
and that the dimensionless shape functions 
are 0(1) in a region where their arguments 
are of order Aqcd • 

• The weak annihilation error (see below) is 
taken as a fixed percentage of the total rate. 

The HFAG average for ICHEP 2006 [21] gives 
the value |V u6 | = (4.49 ± 0.19 ± 0.27) • 10" 3 us- 
ing the BLNP approach, where the first error is 
from experiment and the second from theory. The 
non experimental errors associated with this re- 
sult arc: 4.2% HQ error (for the leading order 
shape function), 3.8% for the combined pertur- 



bative and subleading shape function error and 
1.9% for weak annihilation. 

BLL Approach: The BLL (Bauer-Ligeti- 
Lukc) approach [22] is based on the fact that if 
a low q 2 cut is imposed to eliminate charm back- 
ground, an OPE expansion for the partial rate 
can be constructed, which is suppressed by in- 
verse powers of m c instead of nib- (This is possi- 
ble since for such a cut P+ ~ P_ <~ m c ). In order 
to optimize both the efficiency and the theoretical 
uncertainty, the authors of [22] suggested to use 
a combined Mx — q 2 cut. The main ingredients 
of the BLL approach arc: 

• The OPE is assumed to be valid for the 
combined cut. 

• The LO shape function sensitivity is esti- 
mated by convoluting the tree-level decay 
rate with the difference between the "tree- 
level" shape function model and a delta 
function model for the shape function. 

• The subleading shape function contribution 
is assumed to be small and is not assessed. 

The HFAG average for ICHEP 2006 [21] gives the 
value \V ub \ = (5.02 ± 0.26 ± 0.37) • 10~ 3 using the 
BLL approach, where the first error is from ex- 
periment and the second from theory. The theory 
error includes a 3% contribution from shape func- 
tion sensitivity. 

DGE Approach: The Dressed Gluon Expo- 
nentiation (DGE) approach is advocated by An- 
dersen and Gardi [23]. Similar to the BLNP ap- 
proach it can be applied to various experimen- 
tal cuts, but conceptually it is different from the 
other approaches discussed in this talk. For a 
less technical review of the DGE approach we re- 
fer the reader to [24]. Here we only mention the 
following features of this approach: 

• The decay spectra is approximated by the 
resummed on-shell 6-quark decay spectrum, 
for which the only input parameters are m& 
and a s . 

• Non-perturbative effects associated with 
the meson structure are estimated using 
renormalon analysis. 
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The HFAG average for ICHEP 2006 [21] gives 
the value \V ub \ = (4.46 ± 0.20 ± 0.20) • 10~ 3 using 
the DGE approach, where the first error is from 
experiment and the second from theory. 

Discussion: Let us compare the various ap- 
proaches. As might be expected, the main dif- 
ference between them is the estimate of non per- 
turbative effects and power corrections, beginning 
with the leading order shape function. 

In the BLNP approach the leading order shape 
function is to be extracted from experiment. In 
the current experimental implementation the full 
shape of the photon spectrum is not used and 
only the first two moments of the shape func- 
tion are used to constrain its form. One might 
argue that the possible variations of the photon 
spectrum are already included in the range of the 
HQ parameters, but a verification of this assump- 
tion would be desirable. In the BLL approach 
the shape function sensitivity was included as de- 
scribed above. We should note that the BLNP 
analysis did not find reduced shape function sen- 
sitivity for a combined Mx — q 2 cut. Consid- 
ering the theoretical advances in the control over 
the leading order shape function a reevaluation of 
this issue in the BLL approach is in order. In the 
DGE approach there is no error associated with 
the leading order shape function. Instead there 
is an error from the value of mj, a s and the pa- 
rameter C defined in [23] . (Recently the effects of 
the non perturbative parameter / pv on B — > X s j 
were also considered [25]). It is unclear how this 
C parameter is related to the shape function un- 
certainty in all the other approaches, or to the 
shape function moments. It is also unclear how 
in the DGE approach the OPE result is recovered 
beyond leading order in 1/rnf,, e.g. the terms pro- 
portional to ^ 2 and fifj, when one integrates over 
the kinematical variables. 

Apart from the leading order shape function, 
we have the issue of sublcading shape functions. 
In the BLNP approach the sublcading shape were 
modeled using their moments. In the BLL ap- 
proach the sublcading shape functions contribu- 
tion was assumed to be small and was not as- 
sessed. In the DGE approach there is no er- 
ror from subleading shape functions. Again, it 
is unclear how the parameters C and / PV in 



the DGE approach are related to the sublcading 
shape function error in the other approaches. 

Despite all of these differences, the inclusive 
measurement are all consistent with each other, 
even if we assume that the error bars are under- 
estimated. We have no good explanation for this 
fact. 



5. IMPROVED Kb 

The extraction of \V u b\, impressive as it is, can 
be further improved. We separate the discus- 
sion of improvements that can be implemented 
today, using the currently available theoretical 
tools, and future feasible improvements. 

Improved V u b - Today: In the former class 
we have the treatment of weak annihilation (WA) 
and the weight function approach. 

Weak annihilation appears at order l/mf in 
the OPE. It arises from flavor specific four quark 
operators of the form bTuu Tb, and effects neu- 
tral and charged _B's differently [26]. Currently 
there are only estimates of its magnitude [27]. 
For example a CLEO analysis finds the limit 
T WA /T b ^ u < 7.4% at the 90% confidence level 
[28]. Apart from estimating the WA as a fixed 
percentage of the rate, a different strategy was 
suggested by Lange, Neubert and Paz [5]. Since 
this conference took place in Oxford UK, we can 
follow the Queen of Hearts and summarize this 
strategy as: "Off with its head!" [29]. More con- 
cretely, since weak annihilation is concentrated in 
the region of q 2 = m 2 , if we cut on high q 2 < ^ ax 
(e.g. <7^ ax = (Mb — Md) 2 ) combined with a 
Mx or P + cut to eliminate charm background, 
we would eliminate the WA error. With such a 
cut one loses efficiency but this might be compen- 
sated by the elimination of the WA error. Prelim- 
inary studies in [5] showed that this is indeed the 
case. This method is still waiting for experimental 
implementation. 

The weight function idea is to directly relate 
the photon spectrum of B — > X s j to B — > X u l~v 
spectra without a need to extract the leading or- 
der shape function, and was first suggested in 
[30]. The theoretical input in this approach is 
the weight function, which at leading order can 
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be calculated in perturbation theory: 

w H u ■ J[m b y{P+ - <2Q] ® S(u>) + ... 

In this relation the leading order shape function 
cancels in the ratio. Beyond leading order this 
is no longer the case. BaBar has used a cal- 
culation of a weight function for the Mx spec- 
trum by Leibovich, Low, and Rothstein (LLR) 
[31] to measure \V u b\. The result is \V u b\ = 
(4.43 ± 0.45 ± 0.29) TO" 3 [32], where the first error 
is from experiment and the second from theory. 

A more recent theoretical calculation was per- 
formed by Lange, Neubert, and Paz [33]. In this 
approach a weight function that relates the nor- 
malized photon spectrum in B — > X s j to the P + 
in B — + X u l~v is constructed. The main rea- 
son to use the normalized spectrum is the better 
perturbative convergence of the weight function. 
This weight function contains two loop correc- 
tions from the ratios of jet functions, subleading 
shape function corrections and the known a s /m,b 
corrections. The combined theoretical error is at 
a level of 5%. The approach was generalized by 
Lange [34] and weight functions were calculated 
for an arbitrary B — ► X u l~v spectra. These 
weight functions have the potential to give the 
best extraction of | V u b | • This method is still wait- 
ing for experimental implementation. 

Improved V u b - Future: Apart from these 
already available calculations we can construct 
a wish list for feasible theoretical calculations. 
First, a complete subleading shape function anal- 
ysis for B — > X s j is underway [35] and prelimi- 
nary results were already reported in [12]. With 
the approaching completion of the full order a 2 
corrections to Co for B — > X s -f [36] (see also the 
contribution by T. Hurth at this conference) the 
calculation of H s at 0(a 2 s ) is nearly done [37]. 
In order to construct a full two loop weight func- 
tion a two loop expression for H u is needed which 
is not easy but feasible. The next source of un- 
certainty is the a s corrections for the terms con- 
taining the subleading shape functions. Such a 
calculation is more complicated and connected to 
the calculation of C2 in equation (4) to order a s , 
where the latter would serve as a check to the for- 
mer. Finally it is not completely clear that the 



subleading shape functions cannot be extracted 
from data. With the completion of this wish list 
we would probably reach the boundary of the the- 
oretical accuracy in extracting |V^f,|. 



6. EXCLUSIVE | Kb | AND LESSONS 
FROM LEPTONIC B DECAYS 

\V u b\ can also be extracted from exclusive de- 
cays such as B — * irW (see the contribution by 
L. Gibbons at this conference). In order to do 
so there is a need for a theoretical input about 
the form factor f + (q 2 ). Currently there are two 
different theoretical approaches. The first be- 
ing lattice QCD calculations, valid for q 2 > 16 
GeV, and the second light cone sum rules valid 
for q 2 < 16 GeV. For the first approach HFAG 
[21] cites three sources: an unquenched calcu- 
lation by HPQCD collaboration with \V ub \ = 
(3.93 ± 0.26t^4i) ' 10~ 3 [38] (see also the con- 
tribution by C. Davies at this conference), an 
unquenched calculation by the FNAL collabora- 
tion with |Kb| = (3.51 ± 0.23±g;f ) • 10~ 3 [39], 
and a quenched calculation which can be found 
in [21]. For the second approach there is only a 
single source: a calculation by Ball and Zwicky 
[40] that gives \V ub \ = (3.38 ± 0.12±g;|f) • lO" 3 . 
As can be seen from the above list, the exclusive 
predictions are typically smaller than the inclu- 
sive ones, which was also "historically" the case 
[3,41]. 

Recently the Belle collaboration reported on 
evidence for the (pure) leptonic B — > tD t decay 
[42]. The branching fraction of such a decay de- 
pends on the product of \V u b\ and f B , the B me- 
son decay constant. The product was measured 
to bef B -\V ub \ = (10.1±^(stat)±i; 3 (syst))xl0- 4 
GeV. Using the inclusive value \V ub \ = (4.39 ± 
0.33) • 10~ 3 the decay constant can be extracted: 
f B = . 229±g;glf (stat)±g;g|7 (syst) . This value is 
in good agreement with the unquenched lattice 
value: f B = 0.216 ± 0.022 GeV [43], or the QCD 
sum rule calculations f B = 0.210±0.019 GeV [44] 
and f B = 0.206 ±0.020 GeV [45], supporting the 
value of the inclusive \V ub \- 
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7. CONCLUSIONS 

The last two years have seen an impressive 
improvement in the determination of \V u i,\, es- 
pecially from inclusive decays. This is a result 
of combined experimental and theoretical hard 
work. In this talk we have reviewed some of the 
theoretical work. Further improvement is also 
possible. As we have pointed out, there are the- 
oretical tools that still await experimental imple- 
mentation, namely a cut on high q 2 to eliminate 
weak annihilation, and advanced two loop rela- 
tions between the photon spectrum in B — > X s "/ 
and B — > X u l~~v spectra. Beyond these, more 
theoretical improvement is also feasible. 

The time has now come for a critical com- 
parison of the theoretical approaches to inclu- 
sive | V u b |, namely, comparing the underlying as- 
sumptions, the perturbative and non perturba- 
tive corrections. This is especially important if 
we want to take seriously the 2 a "tension" be- 
tween sin2/3(= sin 20i ) and inclusive \V u b\ 
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